Astronomy & Astrophysics manuscript no. b0952'CII 


©ESO 2013 


March 4,2013 





Resolved [CM] emission in a lensed quasar at z=4.4* 

S. Gallerani 1 - 2 , R. Neri 2 , R. Maiolino 13 , S. Martin 4 , C. De Breuck 4 , F. Walter 5 , P. Caselli 6 , M. Krips 2 , M. 
Meneghetti 7 8 , T. Nagao 9 , J. Wagg 4 , and M. Walmsley 10 ' 11 



(N 

o 

(N 

> 
O 

O 

u 
6 

-t— > 

a 



(N 
> 

in 
m 
o 

in 
o 

(N 



X 



INAF-Osservatorio Astronomico di Roma, via di Frascati 33, 00040 Monte Porzio Catone, Italy 

Institut de RadioAstronomie Millimetrique, 300 rue de la Piscine, Domaine Universitaire, 38406, Saint Martin d'Heres, France 

Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Ave., Cambridge CB3 0HE, UK 

European Southern Observatory, Alonso de Cordova 3107, Vitacura, Casilla 19001, Santiago 19, Chile 

Max-Planck-Institut fur Astronomie, Konigstuhl 17, 691 17 Heidelberg, Germany 

School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK 

INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy 

INFN, Sezione di Bologna, viale Beiti Pichat 6/2, 40127 Bologna, Italy 

The Hakubi Project, Kyoto University, Yoshida-Ushinomiya-cho, Sakyo-ku, Kyoto 606-8302, Japan 
Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy 
Dublin Institute of Advanced Studies, 31 Fitwilliam Place, Dublin 2, Ireland 



Received ; accepted 



ABSTRACT 



We present one of the first resolved maps of the [CII] 158 /im line, a powerful tracer of the star forming inter-stellar medium, at high 
redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-01 15, the host galaxy of a lensed quasar at 
z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two 
[CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over 
a region of about 1 kpc in size near the quasar nucleus, and we infer a star formation surface density > 150 M Q yr _1 kpc -2 , similar 
to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~12 kpc, and located at about 10 kpc 
from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose 
rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These 
observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, 
in the vicinity of the accreting black hole, and in more extended star forming galaxies. 

Key words. Galaxies: high redshift - Galaxies: ISM - quasars: individual: BRI 0952-01 15 - Submillimeter - Infrared: galaxies 



1. Introduction 

The 2 P V2 -> 2 Pi/2 fine structure line of [CII] at 157.74 ^m 
is primarily emitted by gas exposed to ultraviolet radiation 
in photo dissociation regions (PDRs) and is one of the ma- 
jor coolants of the star forming inter-stellar medium (ISM). 
It is the strongest emission line in most galaxies, accounting 
for as much as ~ 0.1 - 1% of their far- i nfrared (FIR) lumi- 
nosity dCrawford et all 119851: IStacev et all 119911: IWright et all 



ll99UlStacev et all 120101: iGracia-Carpio et all 120111) . Given its 
strength, the [CII] line is in principle the most suitable tracer 
to investigate the star forming ISM in galaxies and to iden- 
tify g alaxies th rough the cosmic epochs with (sub-)mm facilities 
(e.g. lMaiolinoll2008l) . 

While in the local Universe the [CII] line is observable only 
from space or airborne telescopes, at high redshift (z > 1) it is 
shifted into the mm/submm windows, and is therefore detectable 
with ground-based observatories. The first detection of the [CII] 
line at high redshift was obtained with the IRAM 30m in SDSS 
Jl 148+5251, one of the most distant qua sars known (z=6.4), fo r 
which the [CII] line is shifted to 1.2 mm dMaiolino et all 12005). 
Since then, an increasing number of [CII] detections at high red- 
shift has been reported, currently for a total of about 20 galaxies, 



* Based on observations made with the IRAM Plateau de Bure 
Interferometer. 



with typical [CII] luminosities L [C in ~ 4.5 x 10 9 - 1.3 x 10 n L o 
(Iono et al., 2006, Maiolino et al., 2009, Hailey-Dunsheath et al., 
2010, Ivison et al., 2010, Wagg et al., 2010, Stacey et al., 2010, 
Cox et al., De Breuck et al., 201 1, Valtchanov et al., 201 1). 

By comparing the observed [CII]/FIR and [CII]/CO lumi- 
nosity ratios with PDR models, it has been possible to con- 
strain the physical conditions of the ISM in star forming re- 
gions of high redshift galaxies and, more specifically, the inten- 
sity of the stellar radiation field, the gas density and the chem- 
ical enrichment (Maiolino et al., 2005, Hailey-Dunsheath et al., 
2010, Ivison et al., 2010, De Breuck et al., 2011, Stacey et 
al., 2010). In the local Universe, the [CII]/FIR luminosity ratio 
drops by an order of magnitude for sources wit h Lhr > 10 11 - 
10 11 5 Lp (luminous infrared galax i es, LIRGs) dMalhotra et all 
12001 iLuhman et all [T9981 I200I iNegishi et all 120011) . This 
trend has been ascribed to various possible effects, such as re- 
duced heating efficiency of the gas in the high UV radiation 
fields, or the presence of a non-PDR contribution to the FIR 
emission, such as might arise from dust-bounded HII regions 
(Kaufman et al., 1999; Luhman et al., 2003; Gracia-Carpio et 
al., 201 1). Mapping the [CII] emission at high redshift can pro- 
vide key information on the spatial disnibution of star forma- 
tion and on the dynamics of primeval systems. However, up till 
the present, it has only been possible to derive estimates of the 
size of the [CII] emitting area towards the host galaxy of the 
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Fig. 1. Map of BRI 0952-01 15 obtained with the PdBI. The synthesized beam of 1.08" x 0.66" is shown in the bottom-left insets. 
Left panel: map of the [CII] line emission region integrated over a velocity range of 300 km s _1 , i.e. -210 < v < 90 [km s -1 ]. 
Contour levels are shown in steps of 2.5<x, where lex = 0.5 Jy km s _1 beam 1 . The filled three-points star, empty square, and filled 
five-points star represent the positions of the components A, B and C, respectively (Table The dashed ellipse denotes the result 
of the fit to the extended component C. Right panel: map of the continuum emission obtained from the line-free channels of the 
3.6 GHz wide spectrum (i.e. v < -210 km s _1 and v > 90 km s ). Contour levels are shown in steps of 2.5tx, where lcr = 0.5 
mJy beam 1 . Dashed lines show the positive contours of the [CII] emission as in the left panel. In both panels, negative contour 
levels at 2.5cr are shown with dotted lines. 



quasar SDSS Jl 148+5251 at z=6.4 (Walter et al., 2009). In this 
object, observations with the IRAM interferometer (PdBI) have 
marginally resolved the [CII] emission and the underlying con- 
tinuum over an area of about one kpc in size, indicating a star for- 
mation surface density of about 1000 M yr~' kpc -2 , similar to 
that observed in the center of some local ULIRGs (e.g. Arp220), 
but about two orders of magnitude larger in area (Walter et al., 
2009). Most likely, such a dense and powerful starburst is trac- 
ing the formation of a massive spheroid in the early Universe. 
In this work, we exploit the new 277-371 GHz band of the 
PdBI to map the [CII] emission in the host galaxy of the quasar 
BRI 0952-01 15. This object was discovered by McMahon et al. 
(1992), and initially identified as a z — 4.5 quasar pair sepa- 
rated by almost 1". After more detailed investigations, the sys- 
tem has been recognized as a lensed galaxy, magnified by a fore- 
ground elliptical at z = 0.632 (Lehar et al., 2000; Eigenbrod 
et al., 2007). The dete ction of BRI 0952-0115 in CO(5^1) by 
iGuilloteau et alj d!999h provided an accurate determination of 
its redshift (z=4.4337) and an estimate of the molecular hy- 
drogen mass (Mh 2 ~ 2 - 3 x 10 9 M Q ). Recently, strong [CII] 
emission has been discovered towards this object using APEX 
(Maiolino et al., 2009). With a peak intensity of 150+30 mJy 
this is one of the brightest [CII] lines detected at z > 4 so 
far. Millimeter and sub-millimeter observations indicate an ap- 
parent far-IR luminos i ty L FT R ~ 1O 13 L for BRI 0952-0115 
dPriddev & McMahonll200lE . However, the [CII]/FIR luminos- 
ity ratio (log(L[di]/LFiR) ~ -2.9, Maiolino et al. 2009) is higher 
than the value observed in local LIRGs and ULIRGs with simi- 
lar infrared luminosities. Additional [CII] observations at high- 
z have confirmed the existence of a large population of distant 
galaxies that, although having Lfir > 10 12 L , are character- 
ized by high [CII]/FIR ratios (Stacey et al., 2010; Ivison et al., 
2010; De Breuck et al., 2011, Valtchanov et al., 2011). These 



high [CII]/FIR - high Lfir - high redshift galaxies are most 
likely extended, scaled up versions of local "normal" starbursts. 
However, it is possible that lower ISM metallicities also con- 
tribute to the [CII] enhancement in high redshift systems as com- 
pared to local systems with high FIR luminosity (De Breuck et 
al.,2011). 



2. Observations 

BRI 0952-0115 [(J2000) RA=09:55:00.1, DEC=-01:30:07.1] 
was observed with the IRAM six-element interferometer (PdBI) 
in 3 observing runs (Jan 15, Jan 16, Jan 17, 201 1) in the extended 
B configuration, and in 2 observing runs (Mar 11, Mar 20, 201 1) 
in the compact C configuration. We tuned the receivers to 349.77 
GHz, which is t he frequency of the [CI I] 157.741 /mi line at red- 
shift z=4.4337 (IGuilloteau et al.Lll999h . 

The beam size resulting from the B+C configurations using 
natural weighting is 1.08" x 0.66", with a position angle of 
20 degree. The lcr sensitivity achieved by the observations is 
0.5 Jy km s _1 beam 1 in the 350 MHz (i.e. 300 km s" 1 ) chan- 
nel associated with the velocity-integrated [CII] emission map 
(see Fig. Q] and Fig. |2), and 0.5 mJy beam 1 in the 3.2 GHz 
(i.e. 3300 km s ) line-free channels used to determine the con- 
tinuum. We also obtain a higher resolution map, characterized 
by a synthesized circular beam of 0.5" and lcr sensitivity of 
0.8 Jy km s _1 beam -1 , by weighting the visibilities with an in- 
verted gaussian taper. 

The spectral setup provides a velocity coverage of ~ 
3600 km s _1 , i.e. ~ 3.6 GHz. A total time of 7.1 hours were 
spent on-source. The weather conditions were acceptable with 
precipitable water vapor 1.5 < pwv < 6.0 mm, which corre- 
spond to a typical atmospheric transmission ~ 0.4 - 0.8. We 
adopt as flux calibrators the following sources: 3C84, 3C273, 
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Fig. 2. [CII] spectrum of BRI 0952-01 15 obtained with the PdBI by integrating over all components (cyan shaded region) and over 
the A+B components (magenta shaded region). The spectrum is rebinned to a resolution of 60 km s _1 (70 MHz). 

Table 1. Properties of the [CII] line observed toward BRI 0952-01 15 with Pd BI (this work) compared with the results obtained with 
APEX (Maiolino et al., 2009) and with the CO(5^1) line observed with PdBI dGuilloteau et allll999h . Uncertainties on the velocity 
integrated fluxes represent statistical errors; for the [CII] lines we report in parenthesis the total errors, which take into account also 
the absolute calibration uncertainties. 



Line 




Instrument 




Vobs 




FWHM 


I 








[GHz] 


[GHz] 




[km s" 1 ] 


[Jy km s~'] 


[CII] CP3/2 


- l P1/2) 


PdBI 


1900.54 


349.776 


4.4336±0.0001 


175±11 


20.1 ± 1.3 (4.2) 


[CII] ( 2 P 3/2 


- 2 P1/2) 


APEX 


1900.54 


349.776 


4.4336±0.0003 


193±32 


33.6±4.9 (7.0) 


CO (5-4) 


PdBI 


576.2679 


106.055 


4.4337+0.0006 


230±30 


0.91+0.11 



MWC349. The absolute flux calibration uncertainty is 20%. The 
data were analyzed by using CLIC and MAPPING (within the 
GILDAS-IRAM packaged). 



3. Results 

The PdBI map reveals a surprisingly complex structure (Fig. [T] 
left panel). This map has been obtained by integrating the line 
over a velocity range of 300 km s , i.e. covering the bulk of 
the line emission at -210 < v < 90 [km s -1 ]. We note that 
the line is slightly skewed towards negative velocities. The [CII] 
map indicates that BRI 0952-0115 is constituted by a compact 
emitting region (labelled A+B in the figure) and a second more 
extended component (labelled C), located about 2" South- West 
of the A+B region. 

The continuum map (Fig. [1] right panel) has been obtained 
by integrating the line-free channels (i.e. v < -210 km s _1 and 
v > 90 km s _1 ). The comparison of the two maps shows that 
the continuum emission and the [CII] line have the same peak 
positions, but display a different morphology. 

Fig. |2] shows the resulting spectrum, rebinned to a spectral 
resolution of 60 kms -1 (70 MHz), and obtained by integrat- 
ing over all components (cyan shaded region). The [CII] line is 
shown on top of a continuum of F cont = 1 1 .4 ± 1 .4 mJy. The line 
was fitted with a single gaussian, centered at z = 4.4336+0.0001 



1 http://www.iram.fr/IRAMFR/GILDAS 



(v = 349.8 GHz), having a FWHM=175 ± 11 km s -1 , and a 
peak intensity F 1 ™"^ = 106.9 + 6.1 mJy. The [CII] line center 
is consistent within lcr w ith that of the CO(5-4) line detected 
bv lGuilloteau et all d 19991) and with the previous APEX obser- 
vations by Maiolino et al. (2009). The [CII] width is also consis- 
tent, within lcr, with the [CII] results by Maiolino et al. (2009), 
while the CO(5^1) line is 25% wider than observed by us, but 
still consistent within 2cr. The PdBI flux is 40% lower than the 
APEX single dish flux, although the two are consistent within 
1.7cr once the absolute calibration errors are taken into account. 
However, if the discrepancy is confirmed with higher accuracy, 
this may suggest that we are missing part of the line flux dis- 
tributed on scales larger than 1". The latter explanation is also 
supported by the fact that the PdBI continuum is ~20% smaller 
than the one obtained through SCUBA observations (McMahon 
et al. 1999), although consistent at ~ lcr. All the results are sum- 
marized in Table Q] 

In Fig. |2 we also plot the bulk of the line emission obtained 
by integrating over the A+B components only (magenta shaded 
region). In this case the line is fitted with a gaussian having a 
FWHM=178 ± 25 km s -1 , consistent with the one obtained from 
the global spectrum, and a peak intensity F ^ /7] = 47.6+5.9 mJy. 

3.1. The lensed quasar host galaxy (component A+B) 

The structure A+B is most likely tracing star formation in the 
vicinity of the quasar nucleus. Moreover, the strongly elongated 
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Table 2. Properties of the BRI 0952-01 15 components compared with optical observations by Lehar et al. (2000): Right Ascension; 
Declination; Flux of each component not corrected for the lensing magnification factor. The uncertainties reported represent statis- 
tical errors. Absolute flux (systematic) errors are not relevant in this context; Flux normalized to component A (our observations); 
Flux normalized to component A (optical observations); Separation from component A (our observations); Separation from com- 
ponent A (optical observations) 





(J2000) RA 


(J2000) Dec 


F [cii] 


(FIFa\cii} 


(F/F A ) opl 


(x — X A )[CII] 


{X %A)opt 




[sec] 


[arcsec] 


[Jy km sec -1 ] 






[arcsec] 


[arcsec] 


A 


09:55:00.125+0.002 


-01:30:06.27±0.04 


5.1+0.5 


1 


1 








B 


09:55:00.088+0.003 


-01:30:06.87±0.05 


3.8±0.6 


0.75 ± 0.20 


0.28 ± 0.01 


0.82 ± 0.06 


0.993 ± 0.003 


C 


09:55:00.021+0.006 


-01:30:08.28±0.10 


9.9+1.3 


2.00 + 0.30 




2.54 ±0.11 





morphology of the structure is due to the blending of the dou- 
ble images lensed by the foreground galaxy. In fact, with the 
synthesized circular beam of 0.5" (see Sec. I3 we can indeed 
resolve the A and B components, as shown in Fig. [3] while com- 
ponent C disappears since its extended emission requires short 
baselines contribution to be detected. For component A we esti- 
mate a radius of < 0.2" (after deconvolving by the beam size), 
i.e. < 1.4 kpc. 

We fit the component A+B in the uv-plane through a circu- 
lar gaussian function with a FWHM = 0.2", finding the coordi- 
nates reported in Table [2] Although the beam is 1.08" x 0.66", 
the signal-to-noise ratio (SNR) of this component is high enough 
that the analysis in the uv-plane allows us to determine the po- 
sition of the two sources with an accuracy of < 0.1". In fact, 
the accuracy of a source position is given by B JV „ f /,/(2 x SNR), 
where B sy „,i, is the synthesized beam size. 

From the high resolution map of Fig. [3]we note that the posi- 
tion of components A and B resulting from the fit in the uv-plane 
with all the baselines (filled three-points star and empty square, 
respectively) do not coincide with the peak of the [CII] emis- 
sion. This may be due to the fact that the diffuse component of 
the [CII] emission has a centroid which is displaced with respect 
to the peak of the [CII] emission. 

We obtain a separation of the lensed components A and B 
smaller than obtained through optical observations by Lehar et 
al. (2000). In our observations, the component A appears to 
be separated by 0.82" + 0.06" from component B, instead of 
0.993" + 0.003", as measured by Lehar et al. (2000). Although 
small, the separation difference is significant relative to the un- 
certainties (Table|2]i. 

We compute the flux associated to the A+B and C compo- 
nents by integrating over regions which follow the 2. 5cr contours 
on the map (see Fig.[TJ left panel), covering ~ 1.7 arcsec 2 and 
~ 2.2 arcsec 2 , in the case of the A+B and C component, respec- 
tively. In particular, the relative flux of B with respect to A is 
the result of the fit to the A+B component. In Table [2j we report 
the flux associated with each component of the [CII] emission. 
One notes both that the [CII] emission is not coincident with 
the optical emission, and that the relative flux of B with respect 
to A resulting from our observations is significantly higher than 
in the case of the HST images. This is likely due to the fact 
that the [CII] emission (tracing star formation) is certainly more 
extended than the optical emission (coming from the accretion 
disk), hence the effect of differential magnification is reduced. 
Alternatively (or in addition) differential dust extinction from the 
lensing galaxy may also play a role. 
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Fig. 3. Map of the [CII] emission of BRI 0952-0115 obtained 
with an inverted gaussian uv-taper to achieve an angular res- 
olution of 0.5" x 0.5". Contour levels are shown in steps of 
lcr = 0.8 Jy km s _1 beam -1 . Negative contour levels at lcr 
are shown with dotted lines. The filled three-points star, empty 
square, filled five-points star, and the dashed ellipse have the 
same meaning as in the left panel of Fig. Q] 



3.2. Lensing model 

In this Section, we discuss the lensing modeQ required to explain 
the structure of BRI 0952-0115 as revealed by optical observa- 
tions and by the PdBI maps. For the mass distribution, we adopt 
the non-singular isothermal ellipsoid (NIE) model. We use the 
GLENS code by Krips & Neri (2004). Eigenbrod et al. (2007) 
identified the lensing galaxy of BRI 0952-01 15 as an early-type 
galaxy at z = 0.632 ± 0.002. 

The model has 7 free parameters: mass (M g ) of the lensing 
galaxy, its core radius (£., which represents the radius at which 
the mass density falls to half of its central value), ellipticity (e) 
of the projected mass distribution and its position angle (PA), 
source position relative to the lensing galaxy in the source plane 
(xs,ys)', source size in the source plane {Rs). We assume that the 
mass distribution of the source is described in terms of a circular 
gaussian function with a FWHM = 0.07". First, we constrain 
the lensing model through the optical position of components A 



We 



concordance A-cosmology with H = 
70.3 km s^Mpc -1 , Cl A = 0.73 and n ra = 0.27 dKomatsu et all 
l201lh 



assume the 



With these cosmological parameters, at z=4.4336, an angular 
distance of 1" corresponds to a proper distance of 6.85 kpc. 
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Table 3. Parameters of the fiducial lensing model: redshift of the source (zs); redshift of the lens (zl)\ mass of the lensing galaxy 
(M g ); ellipticity (e); core radius (£ e ); position angle (PA); source position in the source plane relative to the lens galaxy (xs,ys)', 
source size in the source plane (Rs)', magnification factor of the source (/i Q pt)- 
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Fig. 4. Lensing model of BRI 0952-0115. Left panel: the lens galaxy (at z=0.632) is located at the center of the image. The 
black contours represent the BRI 0952-01 15 position required to match optical observations. The mass distribution of the source is 
described in terms of a circular gaussian function with FWHM~ 0.07". The red curves around the position of the lens denote the 
tangential caustic, while the blue ellipse shows the radial caustic. Right panel: resulting image configuration and critical lines in 
the lens plane. The source is strongly lensed by the foreground galaxy and mapped onto a set of double images: the black contours 
represent the lensed images of BRI 0952-0115. Red crosses represent the relative position of A and B from the lensing galaxy 
location (filled red square), as measured by Lehar et al. (2000) in the optical images. Contour levels are shown at 20%, 50%, and 
80% of the peak. 



and B relative to the lensing galaxy in the lens plane (see Tab. 3 
in Lehar et al., 2000). 

The results of the lensing model are shown in Fig. |4] where 
the panels show the source and the lens planes on the left and on 
the right, respectively. The tangential (red) and the radial (blue) 
caustic lines are displayed in the left panel. The corresponding 
critical lines are shown with the same color-code in the right 
panel. In the lens plane the two red crosses show the position 
of the optical images relative to the lensing galaxy (central red 
filled square), as measured by Lehar et al. (2000), while the black 
contours represent the lensed images of BRI 0952-01 15, as pre- 
dicted by our model. The full set of model parameters is reported 
in Tabled 

Determining the errors on the lensing parameters is a long 
and demanding process, which goes beyond the scope of this 
paper, since we are not interested in the lens itself. In the con- 
text of this paper, the only parameter associated with the lens- 
ing model relevant for our discussion is the lens magnification. 
Hence, here we focus on the determination of the uncertainty of 
the lensing factor. We first compute a large grid over the lens 
shape parameters, and then, for the lensing models which pro- 
vide a good agreement with the data, we calculate a grid over the 
source position to determine Ax 2 < 1 ranges for the maximum 
and minimum values of the magnification factor. By varying the 



coordinates of the source position of ~ 0.002" (see Table 0J, we 
obtain /i opt = 6.7+0.1 for the component A+B, while the C com- 
ponent is not gravitationally amplified. We note that the derived 
magnification factor does not strongly depend on the assumed 
size of the emitting region; in fact, by varying Rs between one 
tenth of the fiducial value (i.e. 0.007") and the radius upper limit 
obtained in Sec.[3T](i-e. 0.2") , // opt changes only by 4%. 



In Sec. 13.11 we have noted that the distance between the 
[CII] components A and B appears to be smaller than that ob- 
tained by Lehar et al. (2000), as shown in Table 2. In the case 
of mm observations, we do not know the relative position be- 
tween the lens and components A and B. Therefore, we con- 
strain the lensing model through the distance in the lens plane 
between components A and B, their flux ratio, and the position 
angle (P.A. = 40° + 2°). We find that the lensing model which 
best describes the properties of the [CII] emitting regions dif- 
fers from the previous one only in terms of the lensing galaxy 
mass, which in this case is 25% smaller. We emphasize that in 
both lensing models the size of component A in the lens plane 
is smaller than 0.2", in agreement with the radius upper limit 
found in Sec. 13.11 By matching mm observations, we obtain a 
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lower value of the magnification factofl i.e. /j mm = 5.8 + 0.7, 
which is however still consistent within \3cr with // opt . 
Further efforts are required to investigate the properties of the A 
and B components in more details. In a forthcoming paper, we 
will present higher SNR and higher angular resolution observa- 
tions of the CO emission in BRI 0952-0115. For the time be- 
ing, we consider as fiducial the amplification factor constrained 
through optical observations, which is characterized by a smaller 
uncertainty. 

3.3. The companion galaxy (component C) 

The C component is more extended than the component A+B. 
We cannot completely rule out the possibility that the C com- 
ponent may result from a residual sidelobe due to a sub-critical 
coverage of the uv-plane; however, we regard this as very un- 
likely, since region C is detected at a confidence level of ~ 7<x. 

We fit the C component with an elliptical gaussian func- 
tion in the uv-plane, finding a size of the emitting region of 
2.5" x 1.1" (~ 17 x 7 kpc) centered at the coordinates reported 
in Table|2] The filled five-points star and dashed ellipse in Fig.Q] 
(left panel) show the results of the fit, while, as mentioned, com- 
ponent C is not detected in the high resolution map of Fig. 4, 
because this map misses the short baselines sensitive to emis- 
sion more extended than 0.5". 

We note that component C is not detected in HST images, 
probably because of its low surface brightness or high extinc- 
tion, which ma kes it undetectable in t he optical. The published 
CO(5^1) map (Gui lloteau et al.L [1999) cannot help to constrain 
the nature of the C component, since the angular resolution of 
these observations is too low (6" X 5") to disentangle the com- 
ponents A+B and C. VLA observations at higher angular resolu- 
tion (~ 2") suggest the presence of CO(2-l) emission extended 
over about 2.5" South- West with respect to the optical quasar po- 
sitions (Carilli et al. in preparation). However, the SNR of these 
maps is not sufficiently high to obtain any definitive conclusion. 

The [CII] line channel maps (Fig. |5j show a complex and 
irregular velocity field, not easy to interpret in terms of simple 
galaxy rotation. Most likely, component C is a companion disk 
galaxy whose velocity field is distorted by the interaction with 
the quasar host galaxy. 

3.4. SFR and E S fr 

We estimate Lfir by using the galaxy template by Polletta et 
al. (2007) of a local ULIRG/QSO (Mrk231) normalized to the 
continuum observed by us at A resl ~ 158 fim. Note that since we 
probe the FIR continuum close to the peak of the SED, the in- 
ferred Lfir is not strongly dependent on the adopted template 
(i.e. on the average dust temperature). For the A+B compo- 
nent, the measured continuum {F conl = 11.4+1.4 mJy) and 
the magnification factor determined through our fiducial lens- 
ing model (fi opt = 6.7 ±0.1) provide a de-lensed luminosity 
L fir B = (1-5 ± x 10 12 L o . If L FIR is associated with star for- 
mation, as in most quasars, even at high redshift (Lutz et al., 
2008), then this luminosit y corresponds to a star formation rate 
SFR = 270 ± 40 M a yr 1 kennicutdfl998h . 

As far as component C is concerned, since the continuum is 
only barely detected, we conservatively use the 3<x upper limit 

3 A lower magnification factor can also be obtained by assuming an 
M/L +y model, as done by Maiolino et al. (2009). 

4 We define the far-IR luminosity as Lp !R = L(42 - 122) yum (Stacey 
etal. 2010). 



on the FIR continuum (i.e. 4.2 mJy). We obtain Lp IR < 3.8 x 
1O 12 L , which corresponds to SFR < 660 M yr -1 . 

Taking into account the upper limit on the size of component 
A mentioned in Sec. 13. II i.e. < 1.4 kpc, and assuming that com- 
ponent A contributes half of the total FIR flux, we infer a SFR 
surface density of Sjp^ <: 150 M Q yr~' kpc -2 . Such a high rate 
of star formation per unit area is typical of the active star form- 
ing region in local ULIRGs (Siebenmorgen et al., 2008; Elbaz et 
al.,2011). 

The compactness of the star formation in the host of the 
quasar nucleus is confirmed by the Ljchj/Lfir ratio. This ratio is 
considered to be a tracer of the radiation field density and there- 
fore a proxy for the compactness of star formation (Stacey et al., 
2010). For the A+B component, we obtain (L [C ii] /L FIR ) A+B = 
(5.3 ± 1.1) x 10~ 4 , typical of local compact ULIRGs (Luhman 
et al., 1998; Luhman et al., 2003). Unfortunately, a resolved CO 
map is still not available to isolate the CO emission of compo- 
nent A+B from component C and, therefore, prevents us from 
locating the individual components on the L[ch]/Lfir versus 
Lco/Lfir diagram, which would allow us to better model the 
physical conditions of the star forming regions in these compo- 
nents. 

We emphasize that the calculations of the SFRD and 
[CII]/FIR ratio do not depend on the magnification factor. By as- 
suming two extreme values of the magnification factor (jx = 3.7 
and /j = 7, i.e. the 3<x lower limit of yu mm and the the 3<x up- 
per limit of ;U p t ), the mean de-lensed FIR luminosity gives a star 
formation rate 150 < SFR < 280 M yr -1 . 

In the companion galaxy C, we can only infer an upper limit 
on the surface density of star formation. Based on the upper limit 
on the SFR and the size measured in the previous section, we in- 
fer Zg FR < 5 M yr~' kpc -2 , which is consistent with moderate 
intensity local starburst galaxies (e.g. M82, NGC253, Arp299). 
The spatial extent of the star formation is also supported by the 
L[cii]/Lfir ratio. More specifically, for the companion galaxy we 
obtain (L[chj/Lfir) c > 2 x 10~ 3 , again typical of extended star- 
burst galaxies. 

This indicates that the interaction between two galaxies at 
high-z can boost both compact star formation and black hole 
accretion in one of the galaxies, while enhanching star formation 
over the whole disk of the other. 



4. Conclusions 

We have presented one of the first resolved maps of [CII] 158/mi 
emission at high redshift. By exploiting the new IRAM PdBI 
receivers at high frequency we have observed BRI 0952-0115, 
a lensed quasar at z=4.4, that is characterized by strong [CII] 
emission. The PdBI map reveals a surprisingly complex struc- 
ture. The [CII] emission can be divided in two main components: 
a compact structure (A+B), associated with the two lensed opti- 
cal images of the quasar, and a second more extended component 
located 2" South- West of the quasar. 

The [CII] emission associated with the quasar (component 
A+B) is clearly resolved in two lensed images as the optical 
quasar images. The [CII] emission associated with the quasar 
nucleus is distributed over a region of 1.4 kpc or less. The con- 
tinuum is also detected in this component, with a de-lensed lu- 
minosity of L^+ B = 1.5 x 10 12 L , implying SFR ~ 270 M yr 1 
. The inferred SFR surface density is Esfr ^ 150 M yr -1 kpc -2 . 
Such a high rate of star formation per unit area is similar to that 
observed in local ULIRGs (Siebenmorgen et al., 2008; Elbaz et 
al., 2011). The compactness of the star formation in the host 
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Fig. 5. Maps of the [CII] emission obtained with the PdBI in 60 km s _1 wide channels. The central velocity of each map (in km s _1 ) 
is indicated in the upper left corner in units of km s _1 . Contour levels are shown in steps of 1.5<x, where lcr = 0.2 Jy km s _1 beam -1 . 
Negative contour levels of 1.5cr are represented through dotted lines. Source C reveals a complex velocity pattern, suggestive of 
rotation disturbed by the interaction with the quasar host galaxy (A+B). 



of the quasar BRI 0952-0115 is also supported by the ratio 
(L[cii] /Lfir) a+b ~ 10~ 4 , characteristic of high UV radiation 
fields typically associated with compact ULIRGs. 

We also detect another [CII] emitting object (C), located 
about 2" South-West of the A+B component. The [CII] emis- 
sion of this component is clearly resolved over about 12 kpc. 
This is likely a companion galaxy, located at about 10 kpc from 
the quasar, in the process of merging with the quasar host. The 
velocity pattern of this component is irregular, suggesting a ro- 
tation field distorted by the tidal interaction with the quasar host 
galaxy. The continuum is only barely detected and indicative of 
a SFR density (< 5 M yr~' kpc 2 ) lower than in the quasar 
host galaxy and more typical of extended starbursting disks. 
The extended and diffuse nature of star formation in this com- 
panion galaxy is further supported by the high Ljchj/Lfir ratio 
(> 2 x 10~ 3 ). 

Galaxy mergers may have quite different effects on the in- 
teracting galaxies. In the case of BRI 0952-01 15 the galaxy in- 
teraction has both boosted compact star formation in the central 
region and black hole accretion in one of the two galaxies, while 
in the other galaxy the interaction has increased star formation 
all over the galaxy disk. 

These results highlights the power of [CII] mapping to in- 
vestigate the nature of quasar host galaxies and, more generally, 
star formation in high redshift galaxies. 
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